Spontaneous vortex phases in ferromagnet-superconductor nanocomposites 
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We report the appearance of spontaneous vortices in lead superconducting films with embedded 
magnetic nanoparticles and a temperature induced phase transition between different vortex phases. 
Unlike common vortices in superconductors, the vortex phase appears in the absence of applied 
magnetic field. The vortices nucleate exclusively due to the stray field of the magnetic nanoparticles, 
which serve the dual role of providing the internal field and simultaneously acting as pinning centers. 
Transport measurements reveal dynamical phase transitions that depend on temperature (T) and 
applied field (H) and support the obtained (H-T) phase diagram. 



PACS numbers: 74.25.Dw, 74.25.Fy, 74.81. Bd 

The interplay between of superconductivity (SC) and 
ferromagnetism (FM) has been attracting the attention 
of the scientific community since the discovery of su- 
perconductivity Recently, there has been a resur- 
gence of this interest due to new phenomena: for exam- 
ple, the increase of the critical current J c in hybrid sys- 
tems containing sub-micron ferromagnetic particles on 
top of type II superconductors @, 0| and the coexistence 
of superconductivity with long-range magnetic order in 
magnetic superconductors JiJ. In a SC-FM hybrid, the 
magnets strongly affect the properties of the supercon- 
ductor, leading to a change in the critical temperature 
T c and critical current J c . Also, they give an oppor- 
tunity to observe new phenomena such as domain-wall 
superconductivity [B| and hysteresis pinning effect 0. 

In order to study the characteristics of the novel spon- 
taneous vortex phases that arise from the interaction be- 
tween vortices and embedded magnetic nanoparticles, we 
fabricated a hybrid system that consists of a 100 nm 
lead film (Pb) containing homogeneously distributed sin- 
gle domain cobalt (Co) particles (mean diameter about 
4.5 nm) with randomly oriented magnetization. 

The samples are produced by the so-called inert- 
gas(Ar) aggregation method with an Ar pressure of about 
10 _1 mbar. It is a co-deposition of Pb and well-defined 
Co clusters directly onto a sapphire substrate without 
buffer and capping layer. The Ar is absorbed by a cry- 
opump at the other end of the cluster chamber and only 
well-defined Co clusters can enter the main chamber. The 
substrate is mounted on a eoldfinger of a rotatable 4 He 
cryostat and is cooled to ~ 40 K during the deposition. 
The samples are deposited at low temperature in order 
to get high quality Pb films. The angle between the ma- 
trix and the cluster beams was 45°. Due to the differ- 
ent beam directions, samples with different Co volume 
fraction can be made within one preparation. The de- 



position rates are controlled by three quartz balances in 
order to monitor the deposition rate at different positions 
of the substrate. Ag contacts for transport measurements 
are pre-deposited and connected to a multi-channel au- 
tomatic measurement system. The typical dimensions of 
the sample were 10mmx3mmx lOOnm. After deposition, 
transport properties in both zero and non-zero magnetic 
field were investigated in-situ with a built-in split-coil 
superconducting magnet (B < 1.2T). 

This method has some advantages: first, the size of 
the Co nanoparticles is tunable and its distribution is 
very narrow, with a standard deviation of less than lnm. 
Second, due to the fact that we do not use any buffer 
or capping layer, any measurement is related exclusively 
to the Pb-Co sample. Third, the magnetic moments of 
the Co particles are randomly oriented so that they have 
zero total magnetic moment in the sample. After depo- 
sition, the samples were annealed at 300 K in order to 
decrease the defect density. The sample of Co clusters for 
the microstructure study was deposited on a carbon foil 
which was mounted on a TEM-catcher. A more detailed 
description of the experimental set-up and operating pro- 
cedures for the same equipment and a similar system of 
Ag matrix and Co particles can be found in the litera- 
ture 0. 

The morphology of the Co nanoparticles was inferred 
by transmission electron microscopy (TEM). As dis- 
played in Fig.[T](a), they are very homogeneous in shape 
with a mean diameter of d~4.5 nm, therefore smaller 
than the vortex core size (see below). Fig. \V[b) charts a 
typical magnetic hysteresis loop measured with a SQUID 
at a temperature just above T c , indicating that the Co 
clusters are indeed ferromagnetic. Similar nanocompos- 
ites have been produced previously However, the 
appearance of vortices in these systems in the absence of 
applied magnetic field has never been observed before. 
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FIG. 1: (Color online) (a) TEM image of the Co nanoclusters. 
(b) Hysteresis loop for a sample with 5% Co volume fraction 
at 8K. 



The nanocomposites were mainly used to investigate the 
pinning effect of magnetic nanoparticles in the supercon- 
ducting phase in the presence of an external magnetic 
field. 

Theoretical calculations predict that due to their stray 
field, ferromagnets inside a superconductor can lead to 
flux bundles in the form of spontaneous vortices, antivor- 
tices, loops and even closed loops [HI In contrast 
with other systems in which superconductivity and ferro- 
magnetism coexist 0, [H[ , the ferromagnetic constituents 
in the present case are small single domain particles with 
randomly oriented magnetization. Both the size of the 
particles and the mean distance between them are smaller 
than the coherence length (£) and penetration depth (A). 
The Co clusters have considerable magnetic moments and 
thus a strong magnetic stray field. If the diameter of the 
particle (4.5 nm) is much smaller than A ~ 80 ran, the dif- 
ference of the stray field between a free magnetic particle 
and a particle in a superconductor is less than one per- 
cent. It means the effect of the superconductivity on the 
stray field of Co particles can be ignored in our case. The 
maximum field strength of this stray field for a single do- 
main Co particle is about 1.3 T at the poles (correspond- 
ing to the Co bulk saturation magnetization of about 1.8 
T) and, therefore, much larger than the critical field of 
B c (0) = 0.1 T for pure Pb. This fact together with the 
well-known suppression of superconductivity due to the 
proximity effect caused by non-superconducting metallic 
particles embedded in a superconductor(cg. Cu particles 
in Pb i3|)makcs a favorable scenario for the formation of 
spontaneous vortices inside the superconducting Pb film. 

Transport measurements give a signature of the pres- 
ence of vortices in a SC and thus can be used to char- 
acterize a vortex state. In order to probe the nucleation 
of the spontaneous vortices in our films, we analyze the 
temperature (T) dependence of the resistivity (p) . 

In Fig.[2l the p-T curves of different samples are shown. 
Fig. [5] (a) gives the p-T curves for an as-prepared film 
(deposition temperature: 40K) with a Co concentration 
of 3.7 v%. One can see a sharp transition at a critical 
temperature T c ^5.5 K. Compared to T c ^7.2 K of an 
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FIG. 2: Resistivity (p) versus temperature (T) for (a) an as- 
prepared film of lead containing about 3.7 v% cobalt nanopar- 
ticles at zero applied magnetic field, (b) a reference sample 
of pure lead prepared by the same method and (c) the same 
sample in (a) but annealed at 300 K. Tg is the temperature 
that characterizes the transition to a zero-resistance state. T* 
is a temperature that indicates the kink in the resistivity. The 
current used in all measurements was I = 0.1 mA. 



undoped pure Pb film, deposited at 40 K and annealed at 
300K (Fig.[2](b)), the Co-doped sample has the expected 
reduction T c caused by the two effects discussed above, 
namely, the proximity effect and the formation of sponta- 
neous vortices. If we compare our T c reduction (AT C ~ 
1.7 K or as 0.5 K/v%Co) with that found in granular 
Cu/Pb structures [l2| where the T c reduction is due to 
the proximity effect only (AT C ~ 0.12 K/v%Cu), we can 
conclude that the effect of spontaneous vortex formation 
in our Pb/Co system probably is the dominating effect in 
the T c reduction. If should be mentioned, however, that 
this reduction in T c due to 3.7 v% Co clusters is much 
smaller than that caused by Co atoms in Pb [l3| ■ There 
a decrease of 2 K/v% Co is observed, which is caused by 
pair-breaking due to spin flip scattering at paramagnetic 
impurities. In the case of our Co clusters, having a su- 
perparamagnetic blocking temperature of ~ 25 K, such 
a process can not occur below T c . 

Annealing the sample at 300 K has dramatic effccts(see 
Fig[2](c)): it not only reduces the resistivity in the normal 
state by a factor of ~5, but also changes the supercon- 
ducting transition. Besides a further reduction of T c , a 
strongly broadened and structured transition is observed, 
having zero resistivity at low temperature but an ohmic 
region between 3.1 and 3.6 K which is followed by a small 
kink and a second increase of the of the resistivity below 
the main transition to the normal state at about 4 K. It 
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is important to point out that although the main feature 
of the p-T curves for several studied sample are similar, 
their exact form and the characteristic temperatures de- 
pends on the Co volume fraction 14 1 . 

In the following we will discuss in detail the structure 
in the p-T curve of the annealed sample shown in Fig. 
[2j From the resistivity of the annealed sample we can 
estimate the superconducting coherence length £ of this 
sample. Since the upper critical field H c2 in our case is 
given by a sum of external applied field and the stray 
field of the Co particles, we avoid estimating the coher- 
ence length £ from H C 2- Instead we assume the dirty 
limit for type II SCs and get the mean-free path from 
the normal state resistivity (see Fig. [21(c)) to determine 
the coherence length from the following expression: 
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where £o = 83 nm is the coherence length of pure lead at 
zero temperature and I is the mean-free path. The val- 
ues of I were obtained from the resistivity using the free 
electron model. The sample mean-free path is / = 3.1 
nm leading to an estimated £ ~ 25 nm for temperatures 
close to 3 K. As discussed above, the ratio between £ and 
A indicates that the sample is a type II superconductor. 
The anomalous p-T curve for H = is a first indica- 
tion of the existence of more than one phase of vortices 
or other forms of trapped flux inside the superconductor. 
Close to the critical temperature, we also see a resistance 
anomaly that is absent in the non-annealed samples and 
it is destroyed by magnetic field. It is known that sim- 
ilar excess resistances were reported in superconducting 
Al nanostructures 15 Jla . 1 1 71 ] and inhomogeneous su- 
perconducting films [18j. Some possible origins for the 
phenomena are normal-superconducting (N-S) interfaces 
induced by dynamic phase slip centers [la . 1 1 €31 ] or non- 
homogeneous distribution of critical temperatures inside 
the sample This anomaly was also reported in a 

system composed of a SC doped with magnetic impuri- 
ties [lj|. We do not know the exact reason for the ap- 
pearance of the resistance peak in our sample. However, 
it is of no relevance for the feature in the p-T curve far 
below the onset of the superconductivity. 

This low temperature H = behavior we see in Fig. [2] 
(c) is similar to the one observed in type II SCs , such 
as Y-Ba-Cu-0 containing random point defects [20[ or 
amorphous SCs such as MoSi [2l| under an applied mag- 
netic field. The curves suggest that at very low temper- 
atures the vortices are organized in a solid state phase. 
Due to the random orientation of the magnetic particles 
and presence of loops, we presume that the vortices are 
highly entangled leading to a kind of vortex glass. Above 
3.1 K, the onset of the ohmic behavior corresponds to 
a phase transition with thermally induced vortex move- 
ment. The kink at T* represents a transition between 
liquids in different pinning regimes or the rearrangement 



and breaking of the vortex loops without depinning from 
the pinning centers. However, from these measurements, 
it is difficult to characterize the phases. We could at- 
tribute this resistivity behavior to either a spontaneous 
creation of vortices that connect to the sample bound- 
aries, some other type of trapped vortex state [lj| or 
flux creep [22| . 

The most convincing evidence of the existence of a 
spontaneous vortex solid (SVS) at low temperatures 
comes from the characterization of the vortex dynamics 
by means of isothermal voltage (V)-current (I) measure- 
ments. Fig. 02a) and Fig. [3jb) show a set of (V-I) curves 
for the same sample of Fig. [2] at H = and tempera- 
tures ranging from 2.26 K to 4.2 K. They indicate the 
existence of spontaneous vortices networks or loops in- 
duced solely by the magnetic nanoparticles' stray field. 
For a better understanding of the relation between the 
p-T and (V-I) curves, we focus on the curves below the 
main transition (Fig. [3](b)) : for low temperatures, the 
voltage is zero for low currents, which indicates the ex- 
istence of a vortex solid. Furthermore, it increases expo- 
nentially with increasing current, a characteristic of both 
vortex glass and flux creep (23| . For intermediate tem- 
peratures, the V-I curves have an S-shape, associated to 
vortex movements with two different activation energies. 
Finally, it increases monotonically with current for tem- 
peratures just below T c , as it is expected for a normal 
vortex liquid. 

On a logarithmic scale (Fig. [3]^a)) a simpler picture 
appears, with a clear separation of two classes of (V-I) 
curves: one for low temperatures, below a critical line Tq 
and the other for high temperatures. The (V-I) curves 
show a positive curvature for higher T and an ohmic be- 
havior for low currents while for low T the curves have 
a negative curvature. The onset of ohmic behavior cor- 
responds to a liquid phase and the exponential decrease 
of the voltage for decreasing values of the current is the 
signature of a vortex glass phase [24[ ■ 

In order to gain a better understanding of the dynamic 
response of this vortex structure, we use a scaling theory 
in the critical region of a vortex-glass transition. In the 
scaling regime, the relation between the electrical field E 
and the current density J is given by 



E ^u{z+2-D) = F ± (\t\- { - D - 1 > J), 



(2) 



where t = (T — Tq)/Tg, Tq is the temperature for the 
vortex glass transition and F± are two universal functions 
for t > and t < [23}]. D is the dimensionality of the 
phase transition, z is the dynamical exponent and v is 
the exponent related to the divergence of the coherence 
length (£ ~ \T — Tg\~ u )- Using this scaling relation with 
D = 3, T G = 3.15 K, v = 3/8 and z = 2.73, the set 
of data from Fig. [3]Ja) ranging from 2.2 K to 3.8 K can 
be scaled to the same universal functions as shown in 
Fig. E^c). This universality lends strong support to the 
presence of a second order phase transition. The values 



4 



of the exponents z and v are smaller but comparable to 
the exponents of the vortex glass transition obtained for 
Y-Ba-Cu-0 containing random point defects [25| and for 
amorphous superconductors [2l|. The difference in the 
exponents gives us a transition with smaller correlation 
lengths and relaxation rates. 

In earlier experiments on embedded magnetic nanopar- 
ticles in SCs, an increase of vortex pinning due to 
the presence of the magnetic particles has been ob- 
served [1, 0, [2^|. In general, these pinning centers are 
extremely efficient and can be originated by various phys- 
ical effects @, d, [27], 28 1 . In our case, we have an even 
stronger constraint: the hybrid systems have an intrinsic 
pinning, since the magnetic particles produce the vortices 
and pin them to their original location. In this sense, the 
difference in the scaling exponents should be related with 
the different type of pinning centers created by magnetic 
nanoparticlcs. We believe that an extension of the scaling 
theory addressing the issue of magnetic pinning is needed 
in order to fully understand this new spontaneous vortex 
solid (SVS). 
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FIG. 3: (Color online) (a) Isotherms of voltage (V) versus 
current (I) in a log-log scale with no magnetic field applied 
and T are 2.26, 2.35, 2.44, 2.53, 2.6, 2.72, 2.8, 2.9, 3.0, 3.1, 3.2, 
3.4, 3.6, 3.8, 4.0 and 4.2 K from right to left, (b) Low-voltage 
region of the same isotherms of V-I. (c) Scaling analysis for 
the same data presented in (a) from T varying from 2.44K to 
3.8K. The exponent s is given by s = v(z — 2 + D). In the 
curves different colors indicate different temperatures, (d) 
Effect of positive and negative magnetic field (varying from 
-0.01 T to 0.01T) in the V-I curves at 2.9 K. For details, see 
main text. 

Now we also can understand the observed sharp super- 
conducting transition of the as-prepared sample shown in 
of Fig. [2] (a): due to the low deposition temperature of 40 
K, the non-annealed sample has a large number of defects 
which results: (i) in a high resistivity of the normal state, 
(ii) a reduction of the £ (£ oc I 1 / 2 ) and (iii) in strong vor- 
tex pinning. Due to the later one the spontaneous vortex 



solid will remain up to the transition temperature T c re- 
sulting in a sharp transition at T c . 

For a further investigation of how the nanoparticlcs' 
magnetization modifies the low temperature vortex state, 
we study the (V-I) curves for B up to 0.01 T, which is 
much smaller than the field necessary to align the mag- 
netic particles (Fig. |Tfb)). The magnetic field is applied 
parallel to the sample surface and we use the following 
procedure: we first align the nanoparticlcs at room tem- 
perature with a magnetic field B = 1 T, in order to 
obtain a net stray field of the Co particles (Bc ), i-e. 
Bco 7^ 0. Next, we field-cool the system and turn off 
the field at a temperature above T c but far below the 
blocking temperature of the particles. We then apply 
the small magnetic field in opposite directions at fixed 
temperature below T c . We find an asymmetric behav- 
ior of the (V-I) curves for different polarities of the field 
(Fig- Eld)). This result can be seen as a superposition of 
two different effects. The relation of the total field (By) 
inside the sample, the external field (B ext ) and the field 
of the Co clusters is Bt+ = B ext + Bc but with reversed 
external field the relation becomes Bt- = B ext — Bc - 
This difference gives rise to a very small shift in the criti- 
cal current that depends on the field direction. A similar 
result was observed in other systems of nanomagnets in- 
side a superconductor Q . A novel and more pronounced 
effect occurs at small currents, when opposite fields play 
the role of moving the system toward the phases above 
or below the vortex glass transition. 




FIG. 4: (Color online) H-T phase diagram of the hybrid sys- 
tem containing Co nanoparticles. SVS represents the spon- 
taneous vortex solid. It is separated by critical line Tg from 
a strongly pinned vortex liquid (PVL). The boundary to the 
unpinned vortex liquid (VL) is T* , the characteristic tem- 
perature of the kink feature seen in figure 2. The last critical 
curve is the main superconducting transition at T c to the nor- 
mal state (N). Inset: resistivity (p) versus temperature (T) 
for various applied magnetic fields from T to 0.3 T (right to 
left). 

For applied fields higher than the coercive field of the 
sample, the (V-I) curves do not depend on the polarity of 
the field. Therefore we analyze the (p-T) (inset of Fig2|) 
and (V-I) data in the presence of an external magnetic 
field and construct an H-T phase diagram. As before, the 
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field is in-plane but now it is strong enough to increase 
the net stray field produced by the Co particles and the 
total field Bt ■ As can be seen in the diagram of Fig. 2J 
it suppresses the superconductivity, shifting the transi- 
tions (To, T* and T c ) to lower temperatures. Together 
with the reduction of the critical temperatures, there is 
a shrink in the SVS phase, showing that the intrinsic 
pinning is stronger than the normal magnetic pinning. 
Finally, we can see that the lines Tq and T* are paral- 
lel to each other. It could indicate that the pinned liquid 
state (PVL) is indeed a phase of partial movement, where 
the vortices disentangle with reorientation and break of 
loops [29]. This non-linear effect could also be a con- 
sequence of the inhomogeneous character of this super- 
conducting phase and different nature of pinning centers 
(magnetic and nonmagnetic) . This can constrain the 
vortex movement and also give rise to regions with dif- 
ferent pinning strengths. 

In conclusion, our transport measurements show the 
signature of novel spontaneous vortex phases in super- 
conducting films with embedded magnetic particles. At 
low temperatures, the spontaneous vortices arc organized 
in a disordered solid state (SVS), similar to a glass. For 
increasing temperatures, the system undergoes a second 
order phase transition to a pinned vortex liquid followed 
by a crossover that resembles vortex depinning. The de- 
pendence of the vortex solid state on the polarity of the 
magnetic field demonstrates that ferromagnetic particles 
inside superconductors can be used for specific forms of 
vortex creation and manipulation. 
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